Measurement of geometric phase for mixed states using single photon interferometry 
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Geometric phase may enable inherently fault-tolerant quantum computation. However, due to 
potential decoherence effects, it is important to understand how such phases arise for mixed input 
states. We report the first experiment to measure mixed-state geometric phases in optics, using 
a Mach-Zehnder interferometer, and polarization mixed states that are produced in two different 
ways: decohering pure states with birefringent elements; and producing a nonmaximally entangled 
state of two photons and tracing over one of them, a form of remote state preparation. 
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When a pure quantum state undergoes a cyclic pro- 
gression, besides the dynamical phase which depends on 
the evolution Hamiltonian, it retains memory of its mo- 
tion in the form of a purely geometric phase factor P, 01 ■ 
This pure-state geometric phase has been experimentally 
demonstrated in various systems such as single photon 
interferometry 01, two- photon interferometry and 
NMR Recently, it has been proposed that fault- 
tolerant quantum computation may be performed using 
geometric phases 0, since they are independent of the 
speed of the quantum gate and depend only on the area 
of the path the state takes in Hilbert space. The next 
step is to investigate the resilience of geometric phases to 
decoherence, and for this a well defined notion of a mixed 
state geometric phase is needed. Some properties of ge- 
ometric phases for mixed states, proposed by Sjoqvist 
et al. l7|, have been recently investigated in NMR 'si- 
Here, we report the first experimental study of geometric 
phase for mixed quantum states with single photons. Due 
to the exquisite control achievable with optical qubits, 
we precisely map the behavior of the phase for various 
amounts of mixture, yielding experimental data in very 
good agreement with theoretical predictions. These re- 
sults are particularly encouraging in light of recent work 
on scalable linear optics quantum computation 0. 

In order to measure a geometric phase, the dynamical 
phase has to be eliminated. It can either be canceled, for 
example, using spin-echo technique for spins in magnetic 
fields yil, or one can parallel transport the state vector 
in order to ensure that the dynamical phase is zero at 
all times. The parallel transport condition for a particu- 
lar state vector \^{t)) is = 0, which implies 
that there is no change in phase when evolves to 
Y^{t + dt)), for some infinitesimal change in time t. How- 
ever, even though the state does not acquire a phase lo- 
cally, it can acquire a phase globally after completing a 
cyclic evolution. This global phase is equal to the geo- 
metric phase, and has its origin in the underlying curva- 
ture of the state space. It is therefore resilient to certain 
dynamical perturbations of the evolution, e.g., it is inde- 
pendent of the speed (or acceleration) of evolution. 

Uhlmann was the first to describe mixed-state geo- 



metric phases in a mathematical context where the paral- 
lel transport of a mixed state is defined in a larger state 
space which purifies the mixed state |l2l| . In this ap- 
proach the number of parallel transport conditions for a 
known N x N density matrix is N"^, but the time evo- 
lution operator U of such a density matrix has only N 
free variables. This approach can only be described in 
a larger Hilbert space with the system and an attached 
ancilla evolving together in a parallel manner [13j . 

Sjoqvist et al. defined a mixed-state geometric phase 
where no auxiliary subsystem is needed 0, 0|- This 
phase can be investigated using an interferometer in 
which a mixed state is parallel transported by a unitary 
operator in one arm; the output then interferes with the 
other arm, which has no geometric phase. The parallel 
transport of a mixed state p — J2k=i Pk\^^) i^l is given by 
{k{t)\k{t)) =0, V A;, i.e., each eigenvector of the initial 
density matrix is parallel transported by the unitary op- 
erator. The resulting N conditions uniquely determine 
the unitary operator and ensure the gauge invariance of 
the geometric phase. One consequence of invariance is 
that each eigenvector acquires a geometric phase 7/c , and 
an associated interference visibility Vk ■ The total mixed- 
state geometric phase factor is then obtained as an aver- 
age of the individual phase factors, weighted by pk- 
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The polarization mixed state of a single photon can be 
represented by its density operator, which can be written 
in terms of the Bloch vector r and the Pauli matrices 
f? = {fa:, CTj,, cr^}, as yO = ^(1 -|- T • a). It represents a 
mixture of its two eigenvectors with eigenvalues ^(1 ±r). 
The length of the Bloch vector r gives the measure of 
the purity of the state, from completely mixed (r — 0) to 
pure (r = 1). For photons of purity r, Eq. becomes 



ue'T, = cos (17/2) -ir sin (f7/2). 



(2) 



where $7 is the solid angle enclosed by the trajectory of 
one of the eigenvectors on the Bloch sphere with corre- 
sponding geometric phase r2/2 (the other eigenvector tra- 
verses the same trajectory, but in the opposite direction. 



leading to a geometric phase — From Eq. ^ we ob- 

tain the visibihty and geometric phase, respectively, 



V = y^cos2 (17/2) +r2 sin^ (^2/2), and (3) 
7g = — arctan (r tan (il/2)) . (4) 

Here 7g is measured in an interferometer by plotting the 
output intensity versus an applied dynamical phase shift 
in one interferometer arm. For pure states, the geometric 
phase given by reduces to half the solid angle (1^/2). 

In our experiment, single-photon states are condition- 
ally produced by detecting one member of a photon 
pair produced in spontaneous parametric downconver- 
sion (SPDC) [13 (we also took data using coherent states 
from a diode laser). Specifically, pairs of photons at 670 
nm and the conjugate wavelength 737 nm are produced 
via SPDC by pumping Type-I phase matched BBO with 
an Ar"'" laser at A = 351 nm. By conditioning on de- 
tection of a 737-nm "trigger" photon (with an avalanche 
photodiode after a 5-nm FWHM interference filter at 737 
nm), the quantum state of the conjugate mode is pre- 
pared into an excellent approximation of a single-photon 
Fock state at 670 nm also with wavelength spread 

6X ^ 5 nm. As shown in Fig. ^ the 670-nm photons are 
coupled into a single-mode optical fiber to guarantee a 
single spatial mode input for the subsequent interferom- 
eter. A fiber polarization controller is used to cancel any 
polarization transformations in the fiber. 

The mixedness of the 670-nm photons is set via two 
different methods [l^- The first uses thick birefringcnt 
decoherers that couple the single photon's polarization to 
its arrival time relative to the trigger 0, 0| . Consider 
a horizontally polarized (|-ff)) and a vertically polarized 
{\V)) photon. Assuming the decoherers delay vertically 
polarized photons relative to horizontally polarized pho- 
tons by more than the photon's coherence length (given 
by A^/AA ^ QOfim), upon detection of the trigger pho- 
ton, \H) will in principle be detected before \ V). Tracing 
over the timing information during state detection erases 
coherence between these distinguishable states; this is 
equivalent to irreversible decoherence [isj . 

To guarantee a pure fiducial state for this method 
of generating mixed states, a horizontal polarizer is 
placed after the polarization controller, followed by a 
half-waveplate (HWP), and finally the decoherers (four 
pieces of quartz of ~3 cm total thickness). By rotat- 
ing the HWP, the state can be prepared in an arbi- 
trary superposition cos 9\H) +sm9\V). The light is then 
sent through the decoherers, effectively erasing the off- 
diagonal terms in the density matrix, resulting in purity 
r — I cos 26*1. In our experiment, the eigenstates of the 
net geometric phase operator are circular polarizations; 
therefore, before entering the interferometer, the quan- 
tum state is rotated with a quarter-waveplate (QWP) 
into a mixture of left {\L) = {\H) + i\V))/V^) and right 
{\R) = {\H) - i\V))/V2) circular polarized light. 



Mixed State Preparation 



Geometric Phase Quantification 
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FIG. 1: Mixed state generation and interferometer to measure 
geometric phase. Mixed states are prepared via two methods: 
1) tracing over the polarization of one photon of a nonmaxi- 
mally entangled polarization state and, 2) using an initial pure 
polarization state with birefringent decoherers that couple po- 
larization to photon arrival time (see dashed box) fv^ . In the 
latter case, tracing over this time prepares a mixed state. 
Half-waveplates at 9i and 62 generate geometric phase but do 
not otherwise alter the transmitted polarization state. Two 
crossed waveplates in the lower interferometer arm give no 
geometric phase, but compensate the optical path difference 
between the arms to achieve high visibility. The interferom- 
eter shape minimizes unwanted polarization changes arising 
from non-normal mirror and beamsplitter reflections. 

Our second method to produce mixed-polarization 
single-photon states, a version of remote state prepa- 
ration is to trace over one of the photons of a 
pair initially in a nonmaximally entangled polarization 
state. This state is prepared using two thin BBO crys- 
tals oriented such that pumping with polarization Op 
produces a variable entanglement superposition state 
cos 0p\HH) + sm0p\VV) 



Here, the first position 



polarization label corresponds to the trigger photon (at 
737 nm) while the second corresponds to its partner (at 
670 nm). A polarization-insensitive measurement of the 
trigger photon prepares the partner in the polarization 
mixed state p670nm = cos'^ 9p\H){H\ + sin^ 9p\V){V\, 
with r = I cos 29p \ . peronm is then transported over the 
single-mode fiber (still with the polarization controller so 
the fiber does not alter the state). As before, a QWP is 
used to rotate the photon's polarization state to a mix- 
ture of \R) and \L) before entering the interferometer. 

After any of the above mixed state preparations, 
the photon is sent into a Mach-Zehnder interferometer 
(Fig. ^ . In the upper arm, the Bloch vector r is evolved 
unitarily using two half-waveplates with optic axes at 9i 
and 92, respectively. The evolution can be illustrated 
(Fig. 121) with one of the eigenvectors of the density ma- 
trix, e.g., traveling along two geodesies going from 
\R) to \L) and back. The trajectory encloses a solid angle 
fl — 4(6*1 — 02) 1131 • The other eigenvector takes the same 
path but in the opposite direction, and therefore encloses 
the solid angle —ft. For mixed states, the length of r is 
reduced, but the same solid angle is subtended. The re- 
sulting evolution fulfills the parallel transport conditions 
for mixed states, and the induced geometric phase is ob- 
tained by substituting il/2 = 29i into equations Q and 
A motorized rotation stage is used to set 9i (to 
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FIG. 2: The sohd angle fl enclosed by the cyclic path of 
one eigenvector of the density matrix. The other eigenvec- 
tor traces the same path but in the opposite direction, thus 
enclosing the solid angle — Jl. O can be varied by adjusting 
9i —62, the relative angle between the optic axes of the two 
HWPs in the geometric phase arm of the interferometer. 

within 0.01°) and thus, the geometric phase. 

To measure 7^ and v, we apply a dynamical phase shift 
in the lower interferometer arm and measure the result- 
ing interference pattern both with a geometric phase (for 
several settings of 9i fill) and without {di — 0). The 
dynamical phase shift is produced with a piezoelectric 
transducer (PZT) on the translation stage on which the 
lower path mirror is mounted. By adjusting the voltage 
across the PZT, the length difference (AL) between the 
arms is varied, giving the probability for the photon to 
exit the interferometer to the detector as 

P(Ai) = i(l + .cos(^~7.)). (5) 

Photons are detected using an avalanche photodiode. To 
conditionally prepare a single-photon Pock state with the 
desired bandwidth, we count only coincident detections 
(within a 4.5-ns timing window) with the trigger detec- 
tor. We estimate the probability of two photons being 
present accidentally during a given coincidence window 
is 3 X 10~® for the decoherer method (using a 4-mm thick 
BBO crystal) and 8 x 10~^ for the entanglement method 
(using two 0.6- mm crystals). Thus the "accidental" co- 
incidence rate (e.g., between photons corresponding to 
different pairs, or from detector dark counts) is negligi- 
ble, and has not been subtracted from the data. 

Data is taken by varying the PZT voltage from 30 to 
70 volts, in 5 volt steps, giving slightly more than one 
period of the interference pattern. At each voltage, data 
is accumulated for 2 s (decoherer method) or 6 s (traced- 
over entangled state method). We plot the number of 
coincidences as a function of PZT voltage, and then fit a 
curve to extract the phase and visibility information for 
each HWP setting |22| . To calculate the phase difference 
due to the geometric phase, we relate the data for each 
HWP setting 61 to the reference data with 9i = (see 
inset in Figs. 3(a) and 3(c)) H^. 

The experimental data are plotted in Figs. |31[a)-|2f), 
along with theoretical curves based on the measured pu- 
rity of the photons. To determine the purity, we mea- 
sure the \H) and \V) components of the mixed state be- 
fore the last quarter-waveplate in front of the interfer- 



ometer. Figs. |3Ia)-(b) show the data for the geometric 
phase and the visibility, respectively, for the experiment 
where the single photons are decohered with thick bire- 
fringent quartz. Figs. EIc)-|3Id) show the corresponding 
data when the mixture is due to entanglement to the trig- 
ger photon. Figs.|nie)-|3If) show results from the coherent 
state, indicating that the data clearly fits the theoretical 
prediction and demonstrates that the single-photon geo- 
metric phase survives the correspondence principle classi- 
cal limit The two geometric phase plots, Fig.|2ta) and 
Fig. EJc) are fiipped along the x axis: in the first setup 
the input states possessed larger right-circular polariza- 
tion eigenvalues, while in the second setup, left-circular 
polarization was dominant |24l | . 

Fig. 3's error bars arise from the fitting program's un- 
certainty estimate of the phase and visibility from the 
raw fringes. This error is consistent with the standard 
error obtained from repeating measurements four times 
to calculate the spread in the geometric phase and visibil- 
ity. We quantify how well the data fits the theory using a 
weighted reduced x^-analysis. For the geometric phases 
(visibilities), we obtain average values of 0.98 (1.36) and 
1.14 (0.94) for the decoherer and entangled state prepa- 
rations, respectively, indicating an excellent fit. Also, the 
values of r retrodicted from our fringe data agree with 
our direct measurements of r within uncertainty. 

We report the first measurement of geometric phases 
for single photons prepared in various polarization mixed 
states, created using two different methods. Specifically, 
we report a novel way of creating decohered one-qubit 
states from entangled two-qubit states, a simple version 
of remote state preparation. Both types of mixed states 
give geometric phase and visibility data in very good 
agreement with the theoretical predictions. Given the re- 
cent advances in linear optical quantum computation , 
and continued interest in geometric quantum computa- 
tion 0, our results indicate that we have a good mea- 
sure of the geometric phase for mixed states, which in 
future work will enable the estimation of fault tolerance 
in geometric quantum computation with linear optical el- 
ements. We also anticipate further experiments on non- 
unitarily evolved mixed states |0| and non-Abelian ge- 
ometric phases |2^ . to ultimately realize a universal set 
of geometric gates for quantum computation |2^ . 
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data is shown inset in (a) {6i = 0° and Oi = 15° for r = 0.81) and (c) (Oi = 0° and Oi = 15° and r = 0.57). 



[4] J. Brendel, W. Dultz, and W. Martienssen, Phys. Rev. A 
52, 2551 (1995); D. V. Strekalov and Y. H. Shih, Phys. 
Rev. A 56, 3129 (1997). 
[5] D. Suter, K. T. Mueller, and A. Pines, Phys. Rev. Lett. 

60, 1218 (1988). 
[6] J. A. Jones et al, Nature 403, 869 (2000); P. Zanardi and 
M. Rasetti, Phys. Lett. A 264, 94 (1999); L. M. Duan, J. 
I. Cirac, and P. ZoUer, Science 292, 1695 (2001); S.-L. 
Zhu and Z. D. Wang, Phys. Rev. Lett. 91, 187902 (2003). 
[7] E. Sjoqvist et al, Phys. Rev. Lett. 85, 2845 (2000). 
[8] J. Du et al., Phys. Rev. Lett. 91, 100403 (2003). 
[9] E. Knill, R. Laflamme, and G. Milburn, Nature 409, 46 
(2001); T. B. Pittman, B. C. Jacobs, and J. D. Franson, 
Phys. Rev. Lett. 88, 257902 (2002); T. B. Pittman et al., 
Phys. Rev. A 68, 032316 (2003); J. L. O'Brien et al. 
Nature 426, 264 (2003). 
[10] A. Ekert et al, J. Mod. Opt. 47, 2501 (2000). 
[11] A. Uhlmann, Rep. Math. Phys. 24, 229 (1986). 
[12] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000), pp. 109-111. 
[13] M. Ericsson et al, Phys. Rev. Lett. 91, 090405 (2003). 
[14] C. K. Hong and L. Mandel, Phys. Rev. Lett. 56, 58 
(1986). 

[15] P. G. Kwiat and B.-G. Englert, in Science and Ultimate 
Reality: Quantum Theory, Cosmology and Complexity, 
J. D. Barrow, P. C. W. Davies, and C. L. Harper, Jr., 
Eds., Cambridge Univ. Press (2004). 

[16] A. J. Berglund, Dartmouth College B.A. Thesis, also 
|quaiit- ph/0010001 (2000); N. Peters et al. Quant. Inf. 



Comp. 3, 503 (2003). 

[17] The classical light is similarly mixed, using an unbal- 
anced polarizing interferometer to separate \H) and \V} 
by more than the ~1 m diode laser coherence length. 

[18] C. H. Bennett et al, Phys. Rev. Lett. 87, 077902 (2001). 

[19] P. G. Kwiat et al, Phys. Rev. A 60, R773 (1999). 

[20] We can fix 02 = 0, since Q depends only on 9i —62. 

[21] The adjustable waveplate must not be wedged, or the 
path overlap will be ruined when the waveplate is rotated. 

[22] Because the voltage applied to the PZT creates a nonlin- 
ear increase in the path length, we fit our data to a si- 
nusoid with nonlinear dependence on the argument. The 
PZT nonlinear fit parameters were determined by mea- 
suring dynamical phase fringes in our interferometer. 

[23] To reduce the effect of dynamical phase drift, the refer- 
ence data is taken both before and after each geometric 
phase-generating HWP setting. The data for each spe- 
cific waveplate setting is used to calculate the geometric 
phase relative to both the initial and final reference data; 
the reported result is the average of the two calculations. 

[24] The r = 0.04 initial state was prepared more right circu- 
lar polarized than left, opposite of the initial states for 
the r = 0.38 and r = 0.57 trials. Therefore, for ease of 
comparison with these curves, the r — 0.04 data curve is 
displayed flipped about the horizontal. 

[25] F. Wilczek and A. Zee, Phys. Rev. Lett. 52, 2111 (1984). 

[26] J. Pachos, P. Zanardi, and M. Rasetti, Phys. Rev. A 61, 
010305(R) (1999). 



